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Abstract 


Multidimensional two-particle Bose-Einstein correlation functions of charged 
hadrons are reported for pp collisions at 2.76 and 7 TeV in terms of dif¬ 
ferent components of the pair relative momentum, extending the previous 
one-dimensional (1-D) analyses of CMS. This allows for investigating the ex¬ 
tension of the source accessible to the femtoscopic correlation technique in 
different directions, revealing a more detailed picture of the emitting source 
in these collisions at increasing energies. The measurements are performed 
for different intervals of the pair average transverse momentum, kr, and for 
increasing charged particle multiplicitiy, N c h. Results in 1-D, 2-D and 3-D 
show a decrease of the fit radius parameters with hr, whereas a clear rise with 
N ch is observed in all cases. In addition, the fit radius parameters at both 
energies show close similarity in size and behavior within the same intervals 
of (N c h, kr)- 
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1 Introduction 


Femtoscopic Bose-Einstein correlations, also known as HBT/GGLP effect, were in¬ 
vestigated in Ref. PI ®J by CMS for pp collisions at y/s =0.9 TeV [3 0], 2.36 TeV 
[3] and 7 TeV [4]. Such phenomenon was discovered by Goldhaber, S. Goldhaber, 
W. Lee and A. Pais (GGLP effect) [T], being the analogous in high-energy collisions 
to a similar method proposed by R. Hanbury-Brown and R. Q. Twiss (HBT effect) 
[2| for estimating angular dimensions of stars. A broad investigation was carried 
out in these studies in terms of the invariant relative momenum Qi nv . In Ref. 0, 
similarly to what was previously observed in e + e - collisions [5], an anticorrelation 
was reported. This result required a study of different fitting functions, as had been 
suggested in the case of small systems, such as the r model [6], which considers 
strong correlations between the space-time coordinates and the momentum compo¬ 
nents of the emitted particles, and was found to describe better the overall behavior 
of the correlation functions. 

A natural extension of those analyses is to investigate the GGLP/HBT correla¬ 
tions with respect to different components of the pair relative momentum, which 
allows for exploring the source sizes in different directions. Such analyses have also 
been studies by other experiments at RHIC and LHC [7, [8j 9-, [10] [IT] [12i IT3] . 
Therefore, the HBT/GGLP correlation is measured in two-dimensions (2-D) as a 
function of the relative momenta along and transverse to the beam direction, qL 
and qr , respectively. In three-dimensions (3-D) the Bose-Einstein correlations are 
studied in terms of qL , qo and qs, these last two obtained, respectively, by pro¬ 
jecting the transverse component q~r in orthogonal directions, i.e., parallel to the 
average transverse momentum of the pair (&t), and orthogonal to both qL and qo- 

Using the same framework as in Ref. mm, the HBT/GGLP effect is further 
scrutinised here for charged hadrons produced in minimum bias events in pp col¬ 
lisions at = 2.76 and 7 TeV with the data collected by CMS at the CERN 
LHC. New 1-D results at both these energies are discussed and then the analysis is 
extended to 2-D and 3-D cases. 

2 Bose-Einstein Correlation measurements 

2.1 Event and Track Selections 

The data sets used in this analysis correspond to minimum bias samples obtained 
in pp collisions at 2.76 and 7 TeV recorded with the CMS detector at the LHC 
[15J. A detailed description of the CMS detector can be found in Ref. [14]. The 
minimum bias sample at yfs = 2.76 TeV was triggered on-line by requiring at least 
one track with pt > 0.4 GeV to be found in the pixel tracker with |? 7 | < 2.4 for a 
pp bunch crossing. Besides, in the offline analysis hadronic collisions were selected 
by requiring a coincidence of at least one Forward Hadronic (HF) calorimeter tower 
with more than 3 GeV of total energy in each of the HF detectors. In the case of pp 
collisions at 7 TeV events were selected by a trigger signal in each side of the Beam 
Scintillation Counter (BSC), coincident with a signal from either of the two detectors 
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indicating the presence of at least one proton bunch crossing the interaction point 
(IP). Collision events were then selected offline by requiring a Beam Pickup for 
Timing for the experiments (BPTX) signal from both beams passing the IP. 

The set of pp collision events at ^/s = 2.76 TeV used in this analysis comprises 
the data collected by CMS in 2013 at the CERN LHC (3.4 million events). At 
yfs = 7 TeV, a combined sample from pp collisions was considered, which uses data 
from three periods of the CMS data taking, i.e., commissioning run (23 million 
events), as well as from the runs 2010A (16 million events) and run 2010B (4 
million events), where the first is almost pileup free, while the later has a non- 
negligible fraction of events with multiple interactions. A filter was used for reducing 
the contamination in case of multiple vertices (the reconstructed vertex with the 
largest number of associated tracks is selected at 7 TeV, while in pp collisions at 
2.76 TeV, an additional primary vertex might be identified as originating from a 
second pp collision by looking at its properties). To assess the related systematic 
uncertainty, an alternative event selection for reducing pileup contamination was 
also investigated by considering only single reconstructed vertex events. 

In the case of pp collisions at 2.76 TeV, three Monte Carlo samples were used. 
For obtaining the BEC results, minimum bias events simulated with Pythia 6 Z2 
Tune CEma were employed, whereas Pythia D6T and Pythia Z2star mm were 
used for estimating the systematic uncertainties related to the choice of Monte Carlo 
tune. Each of them contained 2 million events. For the analysis at 7 TeV about 33 
million Monte Carlo events were simulated using Pythia 6 Z2 Tune. 

The track selection employed in both cases above follows the same criteria as 
in Ref. mm, and are discussed in details in these references, as well as in m- 

2.2 The Bose-Einstein Correlation 

The procedure adopted is the same as described in Refs. 001- Although no par¬ 
ticle identification is considered in this analysis, the contamination from non-pions 
is not expected to be sizeable, as discussed in [3 : . 4J, since pions are the dominant 
type of hadrons in the sample. For each event, the signal containing the Bose- 
Einstein correlations is identified by pairing same charge tracks from the same event 
and distributing them in bins of the relative momentum of the pair, for instance, 
qi 1 = k^ — k^, being fcf the four-momenta of the individual particles in the pair. The 
background distribution or reference sample is formed similarly, by pairing charged 
particles from different events and within the same y range (where the full pseudora¬ 
pidity interval is divided in three subranges Ay, corresponding to —2.4 < y < —0.8, 
—0.8 < y < 0.8, and 0.8 < y < 2.4), as in Ref. [4]. This mixed event technique 
is referred to as “same track density in At/’. A single ratio is then formed, having 
the signal pair distribution as numerator and the reference sample as denominator, 


relative momentum of the pair is defined as Qi nv = — \/ — (&i — fe) 2 ; 

*^sig is the integral of the signal pair distribution of all the events, whereas A/* re f is 
the equivalent in the reference sample. A double ratio technique is then taken with 


with the appropriate normalization, i.e., R= ( 


A W)^W /dQ inv ) 


K 


sig 


(d n 


ref' 


/dQinv) 


The invariant 
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, where Ry[C ' s the s ' n " 


the data and the Monte Carlo single ratios corresponding to the Qi nv distributions, 
in terms of which the Bose-Einstein Correlation (BEC) effect is investigated (31 [4], 

1Z(Q. ) - _R_ - ( djV sig /dQ inv \ j f djv Mc/dQ- nv 

WmvJ ^ dJV rrf /dQ inv j / \dN MC , nf /dQ inV/ 

gle ratio computed with the simulated events generated without BEC. In each case, 
the reference samples for data and simulation are obtained in the same way. This 
double-ratio procedure has the advantage of considerably reducing the sources of 
bias due to track inefficiency and other detector-related effects, as well as other 
Bose-Einstein correlations. 

The GGLP/HBT method reflects not only the quantum statistics of the pair 
of identical particles, but is also sensitive to the underlying dynamics. In particu¬ 
lar, in the case of charged hadrons, the correlation function may be distorted by 
strong, as well as by Coulomb interactions. For pions, the strong interactions can 
usually be neglected in femtoscopic measurements. As in Ref. m. the depletion 
(enhancement) in the correlation function caused by the Coulomb repulsion (attrac¬ 
tion) of equal (opposite) charge pairs in the case of pions is corrected by weighting 
pair-wise with the inverse Gamow factor [18]. This factor, in case of same charge 

G% S (V W ) = 


and opposite charge, is given by G^ s {r] w ) = e ^“_ x , 


with r] w = o^em^Tr/Qinv- where m n the pion mass and Qi nv the invariant relative 
momentum of the pair. 

For performing the multidimensional analysis, the double ratios are investigated 
in terms of the projections of the relative momentum q = ki — k 2 in two or three 
directions. In the 2-D case, the decomposition is made in gz, ( longitudinal com¬ 
ponent, along the beam direction), and qr ( transverse component). In the 3-D 
case, additional projections are considered in the transverse plane, resulting in go 
(outwards), and qs ( sidewards ), respectively along the average transverse momen¬ 
tum of the pair, kr = (hr ± + &t 2 )/2, and orthogonal to it; qo, qs, and gz, are 
mutually orthogonal. This decomposition of the relative momentum of the pair is 
also known as Bertsch-Pratt variables gm ini. The investigations are carried 
out in center-of-mass (CM), i.e., the LHC laboratory frame, as well as in the Local 
Co-Moving System (LCMS), characterized by the frame in which the longitudinal 
component of the pair average momentum (kL = (kL x + ^l 2 )/2 = (k Zl + k Z2 )/ 2) 
is zero. Details about the 2-D and 3-D relative momentum projections, as well as 
the boost to the LCMS can be found in Ref. m 

The parameterizations used to fit the correlation functions in one- (1-D), two- 
(2-D) and three-dimensions (3-D), respectively in terms of Qi nv , (qL,qr) and 
(qs, qL , qo), are listed below. Throughout this analysis h = c = 1 is adopted. 


^(<3inv) = °[ l + Ae (Q inv rt inv r ] (1 + 5 Q inv ), 


( 1 ) 


7 Z(qLqr) — C 


11 + A exp 


q^-Rr + qlRl + 2 qrqLR?LT 


Jjx(l +aq T +Pq L ), (2) 


R(qs,qL,qo) cjl + Aexp 

(1 + aq s + /3q L + 7 qo)- 


q 2 sR 2 s + qlRl + qoR<) + 2 qoqLR 2 Lo 


a/2 


} 


(3) 
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In the above expressions, A is the intercept parameter (intensity of the correlation 
in the smallest bin of the pair relative momentum), C, 5,a,/3 ,7 are constants. The 
exponent a is the Levy index of stability satisfying the inequality 0 < a < 2 . In 
all the above cases, if treated as a free parameter when fitting the double ratios, 
this exponent usually results into a number between the value characterizing the 
exponential (a = 1 ) and the Gaussian (a = 2 ) functions. More details can be found 
in Ref. [23]. For the sake of clarity, we denote the longitudinal component of the 
relative momentum in the CM frame as and in the LCMS, as g£. In particular, 
in the case a = 1 the exponential term coincides with the Fourier transform of 
the source function p(t,x), characterised by a Lorentzian distribution; the radius 
parameters R[ nY , {Rt,Rl), {Rs, Rl, Ro), correspond to the lengths of homo¬ 
geneity fitted to the correlation function in 1-D, 2-D and 3-D, respectively. The 
additional polinomial terms are introduced for accommodating possible deviations 
of the baseline from unity at large values of these variables (long-range correlations), 
as well as for allowing a better quality fit. In the 2-D and 3-D cases, a = 1 leads 
to the so-called stretched exponential function. In Eq. ([ 2 ]), Rt = Rt + t/3t cos <fi 
and Rl = Rl +t/3l where (3t = pk and Pl = originated in the mass-shell 
constraint (q^k^ = 0 -A q° = ^); </> is the angle between the directions of qx 
and kx and r is the source life-time. In Eq. ^ R 2 S = R 2 S , R% = R 2 L + r 2 /3£, 
Rq = Rq + t 2 (3t and R 2 LO = r 2 /^/^. Both in Eq. ([ 2 ]) and ([ 3 J a frame depen¬ 
dent cross-term, respectively proportional to qrqL and qoqL is present. Flowever, 
when the analysis is performed in the LCMS, this cross-term does not contribute 
for sources symmetric along the longitudinal direction. 

2.3 Systematic uncertainties 

Various sources of systematical uncertainties were considered in this analysis, as 
listed in Table [l] the first four being similar to what was discussed in Ref. SHI- 
Two additional studies were also performed: the effect of separating positive from 
negative charges in the single ratios (since (++), as well as (—) charges are added 
in the signal; in the reference sample, besides these two combinations, the (H—) 
case is also added to the sample). The second study added is the effect of pileup 
events, investigated by comparing the results to the case where only single-vertex 
events are considered. The exponential function in Eq. ([l]) (Levy type with a = 1) 
was adopted for this investigation. The total values of the systematic uncertainties 
are calculated by adding the individual contributions in quadrature. The systematic 
uncertainties estimated in 1-D and summarized in Tableware extended to both the 
2-D and the 3-D cases. 


3 Analysis Results 

3.1 One-dimensional results 

The current analysis extends the previous one for pp at ^ = 7 TeV reported in 
Ref. ^ to full data sample, as well as at yfs = 2.76 TeV with full available statistics. 
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Table 1: Spread with respect to the mean values at y/s = 2.76 TeV and 7 TeV 


Systematical Uncertainties 



2.76 TeV 

7 TeV 

Origin of Systematics 

A 

R inv ( fm ) 

A 

R inv ( fm ) 

Monte Carlo tune 

0.032 

0.160 

0.032 

0.160 

Reference Sample 

0.009 

0.047 

0.051 

0.188 

Coulomb Corrections 

0.016 

0.009 

0.018 

0.020 

Track Cuts 

0.014 

0.119 

0.014 

0.119 

Charge Dependence 

0.006 

0.012 

0.007 

0.006 

Pileup filter 

5.0 e-4 

0.011 

0.001 

0.0025 

Total 

0.040 

0.206 

0.065 

0.275 


The corresponding results for single and double ratios are shown in Fig. [T| The fits 
to the double ratios were produced with the exponential function in Eq. 0. with 
a = 1. 

Figure [2] shows results for the intercept parameter A and the invariant radius 
R[ ny from pp collisions at 2.76 and 7 TeV (full statistics), in terms of N c h and 
kr- It can be seen that the results corresponding to the two energies are very 
similar for the different ( N c h , kr) combinations. The intercept A decreases with 
increasing kr and N c h, whereas R[ nY steadily increases with multiplicity and seem 
to decrease with kr, at least for the two largest N c h bins, showing that the lengths of 
homogeneity accessible to interferometric measurements decrease with the average 
pair momentum, as has been previously observed in several different collision systems 
and energy ranges 0 a 0 Dm m ei m .This behavior is expected in the case of 
emitting sources originated from expanding systems. 

In Fig. | (left), R[ nY is investigated as a function of the charged particle mul¬ 
tiplicity, < N ch > (efficiency and acceptance corrected), where the curves are fits 
proportional to N ^ 3 . The results for pp collisions at 2.76 TeV and at 7 TeV are 
consistent with former studies reported in Ref. [4] at 0.9 and 7 TeV, showing a sim¬ 
ilar increase with N^ 3 . Such results also suggest an approximate scaling property 
of the lengths of homogeneity with increasing collision center-of-mass energy. 

As discussed in [4], an anticorrelation (value below unity) was observed in the 
double ratios for values of Qi nv away from the Bose-Einstein peak, whose depth 
was shown to decrease with increasing N c h, for integrated values in This dip 
structure is also observed in the present analysis and its depth is further investigated. 
More details and discussion on the results are presented in Ref. m An additional 
function, 7 Z(q) = C {l + A[cos [(qr^) 2 +tan(mr/4 )(qr a ) a ] e~^ qra ^ a ]} • (1 + Sq), 
was used to fit the data points, which better describes such anticorrelation, as 
discussed in mm- It is based on the so-called r Model [6], which parameterizes the 
time evolution of the source by means of a one-sided asymmetric Levy distribution. 
The dip’s depth [4] is estimated by the difference of the base-line function, C(l+5q), 
and the remaining fit function based on the r model at its minimum, leading to the 
results shown in Fig. [5] (right), where the results for pp collisions at both 0.9 and 7 
TeV from Ref. [4] and the new ones for 2.76 TeV and the full sample at 7 TeV are 
shown together (see m for details). 
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Figure 1: 1-D single ratios as a function of Q[ nY for data and Monte Carlo 

(Pythia 6-Z2 tune) from to pp collisions at 2.76 TeV are shown (left), as well as 
the corresponding double ratio superimposed by the exponential fit (right). 
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Figure 2: The fit parameters A (left) and R[ ny (right) from exponential fits to 
the double ratios are shown in different N c h and kr bins, from pp collisions 
at 2.76 and 7 TeV (full sample). The statistical uncertainties are indicated by 
error bars (in some cases, smaller than the marker’s size), the systematic ones 
by empty (2.76 TeV data) or shaded boxes (7 TeV data). 
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Figure 3: Comparative plots with results in Ref. [4]. Left: Ri ny versus < N c h > 
(acceptance and efficiency corrected), for pp collisions at 2.76 and 7 TeV (fit 
curves are proportional to 7V^ 3 ). The inner error bars represent statistical 
uncertainties and the outer ones, statistical and systematic uncertainties added 
in quadrature. Right: The anticorrelation’s depth, A, versus < N c h >. 
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3.2 Two-Dimensional Results 


The BEC analysis is extended to the 2-D case in terms of the components gz,, Qt of 
the pair relative momentum, with the data samples from pp collisions at yfs = 2.76 
TeV and 7 TeV. It is performed both in the CM frame and in the LCMS, in which 
the cross-term depending on g T gz, in 2-D, or g^go, in 3-D, does not contribute, in 
case of longitudinally symmetric systems. 





Figure 4: The left panel shows double ratios as a function of (g^gr), with 
data from pp collisions at 2.76 TeV (top) and 7 TeV (bottom) in the LCMS 
corresponding to results integrated in all N c h and kr bins. The right panel shows 
the corresponding 1-D projections of the single and double ratios in terms of g£ 
(for \qr \ < 0.05 GeV) and qr (for |g£| < 0.05 GeV). Gaussian, exponential and 
Levy (with a = 1) fit functions are shown superimposed to the data points. 


As an illustration, the double ratios in the LCMS are shown in Fig. [4] as a 
2-D plot (left panel) in terms of (g^gr), and the corresponding 1-D projections 
(right panel) for pp collisions at 2.76 TeV (top), and for 7 TeV (bottom). The 
1-D projections, when plotted in terms of g^, considers only the first bin in q T 
(i.e., qr < 0.05 GeV), and vice-versa, with the data superimposed by the Gaussian, 

















































exponential and Levy (with a = 1) fit functions. 
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Figure 5: The fit parameters obtained with the stretched exponential (Levy- 
type with a = 1) function are shown for pp collisions at 2.76 and 7 TeV in 
the CM frame (top) and in the LCMS (bottom), as a function of (kx) and for 
different N c h bins. The statistical uncertainties are indicated by error bars (in 
some cases, smaller than the marker’s size), whereas the systematic ones are 
indicated by empty (at 2.76 TeV) or shaded boxes (at 7 TeV). 


Analogously to the studies performed in 1-D the double ratios in 2-D were also 
investigated in terms of (qx^qr), in three intervals of the pair average momentum, 
kx, and in different N c x bins. The results from the stretched exponential fit 
(Levy-type with a = 1) [23j to the double ratios, performed both in the CM frame 
(top) and in the LCMS (bottom), are compiled in Fig. [5] The behaviour of the 
directional lengths of homogeneity is very similar in both frames, with Rx (R*x) 
and Rx increasing with charged multiplicity, N c x, and decreasing with the average 
transverse momentum, kx, at least in the larger multiplicity bins, a behaviour similar 
to that observed in the 1-D case and expected for expanding sources. Another 
interesting feature of the data than can be observed in Fig. [5] is that R* L (LCMS) 
> Rx (CM) for the same bins of N c x and kx, suggesting an effect of Lorentz boost 
contraction in the longitudinal length of homogeneity in the CM frame. Regarding 
A, as shown in Fig. [5] no significant sensitivity of the intercept is seen as a function 
of kx- However, within each kx range, A slowly decreases with increasing track 
multiplicity in an similar way in both frames. 

Table [2] collects the values of the radius, R* L ,Rx, and of the intercept, A, 
fit parameters, integrated both in N c x and kx, and obtained with the stretched 
exponential fit (i.e., Levy-type with a = 1). From Table [2] it can be seen that, in 
the LCMS, the rest frame of the pair along the longitudinal direction, R* L « AR T / 3, 
suggesting that the source is longitudinally elongated, at both energies. 

In Fig. [6]the 2-D results for the double ratios versus ( q* L ,qx ) in the LCMS are 
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Table 2: 2-D fit parameters for in the LCMS 



2.76 TeV 

7 TeV 

A 

0.830 ± 0.010 (stat.) ± 0.040 (syst.) 

0.700 ± 0.002 (stat.) ± 0.065 (syst.) 

Rt (fm) 

1.498 ±0.013 (stat.) ± 0.206 (syst.) 

1.640 ± 0.003 (stat.) ± 0.206 (syst.) 

R* (fm) 

1.993 ±0.022 (stat.) ± 0.206 (syst.) 

2.173 ± 0.005 (stat.) ± 0.275 (syst.) 


zoomed along the correlation function axis, which cuts the BEC peak above 1.2. 
Figure [6] also shows the 1-D projections in terms of each of these variables, having 
the complementary one within the first bin, i.e, \qi\ < 0.05 GeV. The results are 
shown in four bins of charged particle multiplicity, N c h, which increases from the 
top left panel to the bottom right one. 



Figure 6: Results obtained in the LCMS for the 2-D double ratios with zoomed 
axes, with the BEC peak cut above 1.2, as a function of (g£, qr) for four charged 
multiplicity bins, N c h, increasing from top left to bottom right. The 1-D pro¬ 
jections in q* L (for \q T \ < 0.05 GeV) and qr (for \q* L \ < 0.05 GeV) are shown 
side-by-side to the corresponding 2-D double ratios. 


3.3 Three-Dimensional Results 

The 3-D correlation function in terms of the variables (qs, q* L , qo) can be visu¬ 
alised through 2-D projections in terms of the combinations (gs,^), (^2,^o), and 
(qo,Qs), w ^h the complementary components within \qo\ < 0.05 GeV, \q$\ < 0.05 
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GeV, \q* L \ < 0.05 GeV, respectively, corresponding to the width of the first bins. 
The 1-D projections in the LCMS of the single and double ratios are shown in Fig. [7] 
for pp collisions at 7 TeV. The points represent the data and the curves the expo¬ 
nential, stretched exponential (Levy-type with a = 1) and Gaussian fit functions. 
The fits are performed to the 3-D double ratios and then projected in the directions 
of qs, qo, similarly to the projections of the data points. 
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Figure 7: The top panel shows the 1-D projections of the 3-D single ratios of 
data and Monte Carlo, in terms of the relative momentum components (qs, q £, 
qo ) analyzed in the LCMS, for pp collisions at 7 TeV; the bottom panel shows 
the corresponding double ratios. The Gaussian, the exponential and the Levy 
(with a = 1) fit functions are shown superimposed to the data points. 


The values of the lengths of homogeneity in the Bertsch-Pratt parameterization 
for pp collisions at 7 TeV, obtained with the stretched exponential fit in the LCMS, 
integrating over all N c h and kr ranges, are summarized in Table [3] together with 
the corresponding intercept fit parameter. Comparing the lengths of homogeneity 
in the 3-D case in the LCMS, from Table [3] it is found that 77^ ~ 1.5 Ro fm and 
R* l ~ 1.2 Rs at 7 TeV. Therefore, the source seems to be more elongated along the 
longitudinal direction in the LCMS also in the 3-D case, with the relation among 
the lengths of homogeneity such as R* L > Rs > Ro- 

The fits to the 3-D double ratios were also investigated in three kr bins (inte¬ 
grating over all N c h). The fit parameters were obtained with Gaussian, exponential 
and Levy (with a = 1) fit functions (when treated as a fit parameter a returned 
values close to unity also in the 3-D case, both in the CM and in the LCMS). The 
results are compiled in Fig. [5] for the data from pp collisions at 7 TeV, showing to 
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Table 3: 3-D fit parameters for pp collisions at yfs = 7 TeV in the LCMS 



7 TeV 

A 

0.568 ±0.002 (stat.) ± 0.065 (syst.). 

Ro (fm) 

1.370 ±0.004 (stat.) ± 0.275 (syst.) 

R s (fm) 

1.784 ±0.004 (stat.) ± 0.275 (syst.) 

Rl (M 

2.105 ± 0.005 (stat.) ± 0.275 (syst.) 


depend noticeably on the type of fit used, the radius parameters being consider¬ 
ably larger in the case of the Levy-type function. The dependence on kr, however, 
seems to be similar for the three fit functions. The Rs fit values are among the 
largest (except for the Gaussian fit) and seem insensitive to kr, in both frames. 
Also Rl seems to be insensitive to kr, and is the smallest radius parameter in the 
CM frame. However, it shows opposite behavior in the LCMS, where its decrease 
with increasing kx is more pronounced, also attaining the largest values of the three 
radius parameters, suggesting an effect related to the Lorentz boost in the longi¬ 
tudinal direction, as in the 2-D case. The Ro fit values are slightly smaller in the 
LCMS as compared to the CM frame, and decrease moderately with increasing kr- 
Its dependence on kr is similar in both frames, although it has a slightly steeper 
decrease with increasing kr in the LCMS. 

The fits to the double ratios were also studied in four N c h bins (integrating over 
all &t)- The corresponding results are shown as a function of ( N c h ) (efficiency and 
acceptance corrected) in Fig. [9] obtained both in the CM frame (top) and in the 
LCMS (bottom). A clear behavior can be seen, common to all fit functions and in 
all directions of the relative momentum components: the fit radius parameters Rs, 
Rl and Ro increase with increasing average multiplicity, indicating an increase in 
the lengths of homogeneity with N c h, similar to what was seen in the 1-D and 2-D 
cases. 

The intercept parameter A was also studied as a function of kr and N c h, both 
in the CM frame and in the LCMS. The corresponding results are shown in Fig. [T0| 
The values of (N c h) shown in the plots were corrected for efficiency and acceptance. 
A moderate decrease with increasing kr is observed. As a function of increasing 
N c h, A first decreases and then seems to saturate. 

In the 2-D and 3-D cases, the values of the longitudinal radius fit parameters 
coincide within the experimental uncertainties m. as expected, since both corre¬ 
spond to the length of homogeneity in the beam direction. From Figures [8] and [9] it 
can seen that there is an approximate scaling of Rl ( R * L ) and Rt with N c h, when 
comparing the results at 2.76 TeV and at 7 TeV. 

In Fig. [TT] the 3-D results for the double ratios in the LCMS (integrated in N c h 
and kr) are shown as 2-D projections in terms of pairs of qs,qL,qo (the comple¬ 
mentary one within \qi\ < 0.05 GeV). The plots are zoomed along the correlation 
function axis, cutting values above 1.2. The corresponding 1-D projections along 
variable qs, qL, and qo (other two variables within \qij\ < 0.05 GeV) are also shown. 
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Figure 8: The radius parameters in the CM frame (top) and in the LCMS 
(bottom), for pp collisions at 7 TeV, obtained from fits to the double ratios 
with three different functions, are shown versus (&t), integrated in N c h . The 
statistical uncertainties are indicated by error bars (in some cases, smaller than 
the marker’s size), whereas the systematic ones are indicated by shaded boxes. 
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Figure 9: Results of fits to the double ratios with three different fit functions 
in the CM frame (top) and in the LCMS (bottom) for pp collisions at 7 TeV are 
shown as a function of (N c h) (efficiency and acceptance corrected), integrated 
in &T- The statistical uncertainties are indicated by error bars (in some cases 
smaller than the marker’s size), whereas the systematic ones are indicated by 
the shaded boxes. 13 
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Figure 10: Results corresponding to the intercept parameter A in the CM frame 
(left panel) and in the LCMS (right panel) are shown on top as a function of (kr) 
(integrated in N c h ), and on the bottom as a function of (< N c h >) (integrated in 
kr)- The values of N c h were corrected for efficiency and acceptance. Statistical 
uncertainties are indicated by error bars (in some cases smaller than the marker’s 
size), whereas systematic ones are indicated by shaded boxes. 



Figure 11: The results for the 3-D double ratios obtained in the LCMS 

are shown, the upper panel corresponding to the 2-D projections in (q^qs), 
(qo,qs), an d integrated in N c h and kx, with \qo\ < 0.05 GeV, 

\ql\ < 0.05 GeV and \qs\ < 0.05 GeV, respectively. The bottom panel shows 1-D 
projections of the same data, with the complementary two variables constrained 
to be within the first bin (i.e., \qij\ < 0.05 GeV). 
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4 Summary and conclusions 

The analysis discussed here extends the 1-D results reported in the two BEC publica¬ 
tions in El 01 by measuring the correlation functions in terms of different components 
of the pair relative momentum, as is usually studied by other experiments. This al¬ 
lows to investigate the extension of the source accessible to the correlation technique 
in different directions. Two main projections are considered: in two-dimensions (2- 
D), the femtoscopic correlation is investigated as a function of the variables qL and 
qr, and in three-dimensions (3-D), as function of qL,qs,qo, being qL the same 
as in the 2-D case. For achieving this purpose, minimum bias events produced in 
proton-proton collisions at 2.76 and 7 TeV (full data sample) are scrutinised in de¬ 
tail, as if by means of a magnifying lens. At 7 TeV, the full data sample is used and 
the corresponding results are compared with the ones at the same energy recorded 
during the commissioning run at the LHC in 2010. This analysis also extends the 
measurements of the BEC correlations to the full minimum bias sample from pp 
collisions at 2.76 TeV collected in 2013, which is a very important baseline for the 
measurements of this second order interferometry in PbPb collisions at the same 
energy per nucleon. In 1-D are the results from both energies are compared with 
the ones in Ref. El 01 at lower energies, as well as at yfs = 7 TeV, recorded during 
the commissioning run at the LHC in 2010. In particular, comparisons showed that 
i^iny steadly increases with the charged multiplicity proportionally to nT- 

The measurements were performed both in the collision center-of-mass (CM) 
frame and in the Local Co-Moving System (LCMS), where the average longitudinal 
momentum of the pair is zero. In the 2-D case, for integrated values of N c h and 
kr, the lengths of homogeneity in the LCMS suggest that the source is elongated 
along the beam direction, i.e., R* L > Rt ■ In the 3-D case, it was found that 
R* l > Rs > Ro ■ In addition, it can be observed that the fit values for the 
longitudinal radius parameter, Rl are consistent in 2-D and in 3-D cases, as should 
be expected, since they correspond to the length of homogeneity in the longitudinal 
direction in both cases. This conclusion is attained with respect to the Rl fit 
parameter in the CM frame, as well as R* L in the LCMS m 

The anticorrelation observed in 1-D and reported in Ref. [4], was also observed 
in minimum bias pp collisions at 2.76 TeV and further investigated here with the full 
statistics at 7 TeV. These new 1-D results are also compared with those in Ref. [4], 
verifying a consistent behavior both for the invariant radius parameter and for the 
dip’s depth measurements at 2.76 and 7 TeV. Such comparisons showed that the 
dip’s depth decreases with increasing N c h- 
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